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Abstract 

We consider the bifurcation problem u"+Xu = N{u) with two point boundary 

conditions where N{u) is a general nonlinear term which may also depend on 

the eigenvalue A. We give a variational characterization of the bifurcating 

branch A as a function of the amplitude of the solution. As an application 

we show how it can be used to obtain simple approximate closed formulae for 

the period of large amplitude oscillations. 
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Different physical systems reduce to the consideration of a one dimensional equation of 
the form u" = F{u,fi), where the force F depends on u and one or more parameters fi. In 
these systems the energy E = ^u'"^ + V{u) is a constant of the motion and the qualitative 
dynamics can be understood from the analysis of the phase plane and the stability of the 
fixed points. Among the systems that fall in this class we may cite the buckling of a thin rod, 
the pendulum, steady state solutions of the reaction diffusion equation and numerous others 
IHJ^^,^] (Given the large number of references in this area we refer to some of the standard 
textbooks only). We shall assume that the system has an equilibrium point at m = 0. For 
small deviations from the equilibrium point the behavior of the solution is determined by 
the linear equation u" + Am = with suitable boundary conditions. An important feature of 
the linear problem is that the eigenvalue A does not depend on the amplitude of the motion. 
For large deviations from the equilibrium, the full problem, which we write as 

u" + \u = N{u, A) (1) 

must be considered. A solution to this equation and the analytic determination of the 
eigenvalue A is possible only for special forms of the nonlinearity A^. When an analytic 
solution is not possible, one can resort to obtaining a numerical solution, or use perturbation 
theory around the linear problem. The eigenvalues of the linear part of Eq.(|l]) correspond to 
the points where bifurcation occurs in the nonlinear problem. There exist different methods 
of perturbation theory applicable to these systems An important property of the 

solution in the nonlinear problem is that the eigenvalue A depends on the amplitude of the 
solution. 

The purpose of this article is to show that the eigenvalue derives from a variational prin- 
ciple. The main tool used in the derivation of this result is a new variational characterization 
of the linear problem which, for the sake of clarity, we consider first. 

We wish to solve the linear problem 

u" + Xu = (2a) 

subject to 
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u{0) = 0, = (2b) 

where the boundary conditions we have chosen are the ones appropriate for a wide range of 
problems. We refer to other choices at the end. We will denote u{0) = Um- The positive 
solution to this problem is given by MmCos(7rx/2) corresponding to the lowest eigenvalue 
A = (7r/2)^. There exists a variational principle for this equation, the Rayleigh-Ritz principle 
which provides upper bounds on A. This principle however cannot be extended to calculate 
the eigenvalue of nonlinear problems. We shall construct a new one whose extension to 
nonlinear problems will be straightforward. 

Let g{u) be an arbitrary positive function such that g{0) = and g'{u) > 0. Multiplying 
Eq.(PaD by u'g{u) and integrating in x we obtain, after integrating by parts, 

— - / {uYg\u)u'dx + \ / ug{u)u dx = 
2 Jo Jo 

where the boundary terms vanish since u'{0) = and 5f[u(l)] = g{0) = 0. From the equation 

and the boundary conditions it follows that < in (0, 1) so in the second integral we may 

use u as the independent variable. We obtain then the identity, valid for any admisible g, 



A / ug{u)du = —- I (u) g'{u)u'dx. (3) 



Consider now, for a fixed g the functional 

JM = fwfg\v)v'dx ^vy)dx (4a) 

I Jo 2, Jo 

defined for functions v which satisfy 

v{Q)=Um, ^^(1) = 0, and f ' < in (0,1). (4b) 

Then Jg[v\ > and for fixed g it has a unique minimum attained at v = v. A rigorous proof 
will be given elsewhere. The minimizing function v satisfies the Euler-Lagrange equation 
for Jg. Since $(f,f') does not depend on x the Euler-Lagrange equations can be integrated 
once to yield $ — v'd^/dv' = constant. Since ^{v,v') = —hv^g'{v) we get 



vTg'iv) = -K. (5) 



Then, given g we may determine v by integrating Eq.(|) subject to the boundary conditions 
(^bl). We have then, for fixed 

Jg[v] > min Jg[v] = Jg[v] = ^K{um) (6) 

where the dependence of K on Um is obtained through the boundary conditions on v. 

Before going any further, consider for example g{v) = v. Then, the minimizing v, v, 
satisfies v'{x) = —K^^^. Imposing the boundary conditions we get v = K^^^{1 — x) provided 
that K = u^. For this choice of g, therefore, Jg > 

Now we go back to the problem under consideration. We have from equations (^, (^a]) 
and d^), 

A / ug{u)du = — / {u') g'{u)u'dx 
Jo 2 Jo 

> Jsi^] = -I ['{vTg'ivWdx (7) 
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which gives an upper bound on A. We have shown that for any given g{u) (with g{0) = 
0, g'{u) > 0), 

-2£-ug{u)du 

When will equality hold in Eq.(^? From Eq.(|^) we see that the equality will hold when g 
is chosen in such a way that v coincides with u. For this g which we will call g we will have 

u''g'{u) = -K. (9) 

Since m is a solution of the differential equation (2a) we know that 

+ l^u' = E = hul (10) 

from where it follows m'^ = A-u^(l — [u/umY). Replacing this in Eq.(||) we obtain 
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[{\ul){l-u^/ul)f'- 
4 



This equation can be integrated to yield (omitting overall multiplicative constants) 

= (1 -(«./«.„.)=) v^- 

Our final result for the linear problem is then the following variational characterization for 
the lowest eigenvalue, 

A = max J— ^^^^4:-, (12) 

where the maximum is taken over all positive functions g such that g{Q) — 0, g' > 0. The 
maximum is attained when g — g. 

It is straightforward to evaluate the integrals for g = g and verify that we obtain the 
correct result. Solving v'^g' — —K with v(0) = li^, v{l) = we obtain K — (7r/2)^M^ and 
/ ug{u)du — (l/2)(7r/2)it^ and the exact value for A is obtained. 

We now consider the nonlinear problem 

u" + \u^N{u) (13a) 

subject to 

u'{Q) = 0, u{l) = (13b) 

As before we denote w(0) = Um- Again multiplying by u' g{u), where g is as before, we 
obtain, after integrating 

X ug{u)du^ N{u)g{u)du-- {u'fg'{u)u'dx. (14) 
Jo Jo 2 Jo 

where in the term involving A and the nonlinearity we have used u as the independent 
variable. The second term in the right side is just the linear term considered above so we 
have 

1 Jo"- N{u)g{u)du + K{um) . . 

^-2 !^-ug{u)du 

and, as before the maximum is attained when g is chosen in such a way that v = u. Now 
the equation for g becomes. 



where V{u) is the potential. Here the potential V = |Am^ — Jq N{y) dy and E = |An^ — 
jum ]\f^y"j dy Thls expresslon for g cannot be integrated in general due to the nonlinear terms 
in V{u) that arise from N{u). However, as in the linear case, the maximizing g exists and 
we obtain our main result 

/o"'" N{u)g{u)du + Kju^) 

A — max — — — r — I 

2 J^'-ug{u)du ^ ^ 

where the maximum is taken over all positive functions g such that g{0) = and g'{u) > 0. 
The maximum is attained for g = g. 

Since in general we will not be able to integrate Eq.(|16D the above principle gives lower 
bounds which can be arbitrarily close to the exact value. 

As an example we shall use as a trial function the function g which gives the correct 



eigenvalue for the linear problem g{u) = — [u/umY. With this trial function we 
obtain from ([T^) the lower bound 

/7r\2 4 N{u]u , 

A > 77 + / I ' du (18) 



2^ ■ TTM^^O ^l-u^/ul^ 

which is valid for any N{u). 

As a first example consider the equation 

7 1 

u" + Am = + -u^ 

5 2 

with u{0) = Um, u'{0) = 0, u{l) = 0. The integral in Eq.(]T8|) can be done easily. We obtain 

, ^ /7r\2 4 /2 3 TT 7 2 1 8 3 \ 

A > - + - -Urn ui H ui . 

- \2j TT V3 8 2 5 2 15 

In Fig. 1 we show a plot of nm(A). The solid line corresponds to the exact value 
obtained from the numerical integration of the equation. The dashed line corresponds 
to the lower bound given by the formula above. As it can be seen in the figure, 
the lower bound gives a very accurate value for amplitudes up to Um ~ 2. For 



larger values there is a small discrepancy which can be minimized by a more adequate 
choice for g. The formula given in any case is a close enough lower bound even at 
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For this equation x = is a stable equilibrium point and we search for the period of 
oscillation of perturbations from the equilibrium. This is an exactly solvable equation and 
the period can be given in terms of elliptic functions. When the nonlinearity is small, 
perturbation theory can be used to calculate the period or frequency of oscillations. The 
frequency obtained assuming 5 ^ 1 is 

27r 3 15 

ujpert — — — 1 + -Sa^ — T^^^CL^ + higher order terms 
1 o zoo 

where a is the amplitude of the oscillations. In order to make use of the variational result 
we recall that since the potential has the symmetry x — >■ —x we may consider a quarter 
period of the solution which we chose as the lower right quadrant in phase space {x,x). In 
this region of phase space we have x{0) — 0, x{0) — a, and x{T/4) — 0. Then, introducing 
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the scaled time variable r = 2ujt/TT we obtain 

cPx / ^ \ 

with 

a;(0) = a, x(0) = 0, x(l) = 0. 
We identify A = (7r/2co')^, A^(x) = — 5Ax^ and apply Eg. (P^) . We obtain the bound 

which written in terms of uo gives, 



CU.ar < yi + ^5a^ (19) 

This bound is valid for oscillations of any amplitude. In Fig. 2 we show the results 
for 5 = 0.1. The solid line gives the exact solution, the dashed line the variational 
bound and the dot-dashed line the perturbation result including terms up to order 6"^. 
Again we see the close agreement up to fairly large amplitude of the variational expres- 
sion. As in the previous example a different trial function would give a close estimate at 
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larger amplitudes. 

Fig. 2 
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In this last example we were able to reduce to a quarter period of the solution due to 
the symmetry of the potential. In the absence of such symmetry we should consider half 
the period, that is the lower half of the phase plane. The relevant bifurcation problem to be 
considered in this case is 

u" + \u = N{u) 

subject to 

m(0) = Urn, m'(0) = 0, U\l) = 0. 

One can obtain a similar variational principle in this case too. The details are slightly more 
involved but entirely similar to what we have shown. 

To conclude, we have shown that nonlinear eigenvalue problems for hamiltonian systems 
derive from a variational principle. The main tool in the derivation of this result is a new 
variational formulation of the linear problem. The variational principle may be used to 
calculate the eigenvalues as accurately as desired by suitable choice of trial functions. A 
simple estimate, using as a trial function the optimal trial function for the linear problem, 
gives reasonably accurate results for small to medium amplitude solutions, and it gives a 
good rigorous bound for larger amplitudes. The method that we have used here can be 
applied to other problems. In particular it can be shown to be a reformulation in real space 
of the method used previously by us to obtain the speed of fronts of a reaction-diffusion 
equation in one dimension for arbitrary nonlinearities. Due to the nature of the front problem 
though, its formulation directly in phase space is simpler. It is an open question whether 
this approach can be extended to treat limit cycles of nonlinear oscillators. 

This work has been partially supported by Fondecyt project 1960450. R.B. was supported 
by a Catedra Presidencial. 
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